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Abstract. - We study coupled irreversible processes. For linear or linearized kinetics with mi- 
croreversibility, x = Kx, the kinetic operator K is symmetric in the entropic inner product. This 
form of Onsager's reciprocal relations implies that the shift in time, exp(A't), is also a symmetric 
operator. This generates the reciprocity relations between the kinetic curves. For example, for 
the Master equation, if we start the process from the ith pure state and measure the probability 
Pj{t) of the jth state {j 7^ i), and, similarly, measure Pi{t) for the process, which starts at the jth 
pure state, then the ratio of these two probabilities Pj{t)/pi{t) is constant in time and coincides 
with the ratio of the equilibrium probabilities. We study similar and more general reciprocal 
relations between the kinetic curves. The experimental evidence provided as an example is from 
the reversible water gas shift reaction over iron oxide catalyst. The experimental data are ob- 
tained using Temporal Analysis of Products (TAP) pulse-response studies. These offer excellent 
confirmation within the experimental error. 



Introduction. — 

A bit of history. In 1931, L. Onsager [HIS] gave 
_ the backgrounds and generalizations to the reciprocal re- 
• lations introduced in 19th century by Lord Kelvin and 
H. v. Helmholtz. In his seminal papers, L. Onsager men- 
' tioned also the close connection between these relations 
[ and detailed balancing of elementary processes: at equilib- 
. rium, each elementary transaction should be equilibrated 
I by its inverse transaction. This principle of detailed bal- 
■ ance was known long before for the Boltzmann equation 
I [3]. A. Einstein used this principle for the linear kinetics 
1 of emission and absorption of radiation [4]. In 1901, R. 
Wegscheider published an analysis of detailed balance for 
chemical kinetics [5]. 

The connections between the detailed balancing and 
Onsager's reciprocal relations were clarified in detail by 
N. G. V. Kampen [6j. They were also extended for vari- 
ous types of coordinate transformations which may include 
time derivatives and integration in time [7]. Recently, [8], 
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the reciprocal relations were derived for nonlinear coupled 
transport processes between reservoirs coupled at meso- 
scopic contact points. Now, an elegant geometric frame- 
work is elaborated for Onsager's relations and their gen- 
eralizations ^ . 

Onsager's relations are widely used for extraction of ki- 
netic information about reciprocal processes from experi- 
ments and for the validation of such information (see, for 
example, |10[): one can measure how process A affects 
process B and extract the reciprocal information, how B 
affects A. 

The reciprocal relations were tested experimentally for 
many systems. In 1960, D.G. Miller wrote a remark- 
able review on experimental verification of the Onsager 
reciprocal relations which is often referred to even now 
[llj . Analyzing many different cases of irreversible phe- 
nomena (thermoelectricity, electrokinetics, isothermal dif- 
fusion, etc). Miller found that these reciprocal relations are 
valid. However, regarding the chemical reactions. Miller's 
point was : "The experimental studies of this phenomenon 
... have been inconclusive, and the question is still open 
from an experimental point." 

According to Onsager's work [1], the fluxes in chemi- 
cal kinetics are time derivatives of the concentrations and 
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potentials are expressed through the chemical potentials. 
The fluxes (near equilibrium) are linear functions of poten- 
tials and the reciprocal relations state that the coefficient 
matrix of these functions is symmetric. It is impossible to 
measure these coefficients directly. To find them one has 
to solve the inverse problem of chemical kinetics. This 
problem is often ill-posed. 

Such a difficulty, appearance of ill-posed problems in the 
verification of the reciprocal relations, is typical because 
these relations connect the kinetic coefficients. Sometimes 
it is possible to find them directly in separate experiments 
but if it is impossible then the inverse problem arises with 
all the typical difficulties. 

In our work we, in particular, demonstrate how it is 
possible to verify the reciprocal relations without the dif- 
ferentiation of the empiric kinetic curves and solving the 
inverse problems, and present the experimental results 
which demonstrate these relations for one reaction kinetic 
system. For this purpose, we have to formulate the recip- 
rocal relations directly between the measurable quantities. 

These reciprocal relations between kinetic curves use 
the symmetry of the propagator in the special entropic in- 
ner product. A dual experiment is defined for each ideal 
kinetic experiment. For this dual experiment, both the ini- 
tial data and the observables are different (they exchange 
their positions), but the results of the measurement is es- 
sentially the same function of time. 

The structure of the paper. We start from the classi- 
cal Onsager relations and reformulate them as conditions 
on the kinetic operator K for linear or linearized kinetic 
equations x = Kx. This operator should be symmetric 
in the entropic inner product, whereas the matrix L that 
transforms forces into fluxes should be symmetric in the 
standard inner product, i.e. Lij = Lji. The form of recip- 
rocal relations with special inner product is well known in 
chemical and Boltzmann kinetics [ IWH] . They are usually 
proved directly from the detailed balance conditions. Such 
relations are also universal just as the classical relations 
are. 

Real functions of symmetric operators are also symmet- 
ric. In particular, the propagator exp(ii't) is symmetric. 
Therefore, we can formulate the reciprocal relation be- 
tween kinetic curves. These relations do not include fluxes 
and time derivatives, hence, they are more robust. We 
formulate them as the symmetry relations between the 
observables and initial data (the observables-initial data 
symmetry). 

A particular case of this symmetry for a network of 
monomolecular chemical reactions or for the Master equa- 
tion, which describe systems with detailed balance, seems 
rather unexpected. Let us consider two situations for a 
linear reaction network. 

1. The process starts at the state "everything is in Aq'\ 
and we measure the concentration of Ar. The result 
is c"(t) ("how much Ar is produced from the initial 



2. The process starts at the state "everything is in Ar" , 
and we measure the concentration of Aq. The result 
is Cq{t) ("how much Aq is produced from the initial 
Ar") (the dual experiment). 

The results of the dual experiments are connected by the 
identity 

eq — cq J 
i^r 

where c are concentrations and c°'i are equilibrium con- 
centrations. 

The symmetry with respect to the observables-initial 
data exchange gives the general rule for production of the 
reciprocal relations between kinetic curves. 

Many real processes in chemical engineering and bio- 
chemistry include irreversible reactions, i.e. the reactions 
with a negligible (zero) rate of the reverse reaction. For 
these processes, the micro-reversibility conditions and the 
backgrounds of classical Onsager relations are not appli- 
cable directly. Nevertheless, they may be considered as 
limits of systems with micro-reversibility when some of 
the rate constants for inverse reactions tend to zero. We 
introduce the correspondent weak form of detailed bal- 
ance, formulate the necessary and sufficient algebraic con- 
ditions for this form of detailed balance and formulate the 
observables-initial data symmetry for these systems. 

The experimental evidence of the observables-initial 
data symmetry is presented for the reversible water gas 
shift reaction over iron oxide catalyst. The experimental 
data are obtained using Temporal Analysis of Products 
(TAP) pulse-response studies. These offer excellent con- 
firmation within experimental error. 

Two forms of the reciprocal relations: forces, 
fluxes and entropic inner product. Let us consider 
linear kinetic equations or kinetic equations linearized near 
an equilibrium x^"^ (sometimes, it may be convenient to 
move the origin to x'^'^): 

x^Kx. (1) 

In the original form of Onsager's relations, the vector 
of fluxes J and the vector of thermodynamic forces X are 
connected by a symmetric matrix, J = LX, Lij — Lji. 
The vector X is the gradient of the corresponding thermo- 
dynamic potential: Xi = d^/dxi. For isolated systems, 
<I> is the entropy. For other conditions, other thermody- 
namic potentials are used. For example, for the constant 
volume V and temperature T conditions, $ is —F/T and 
for the constant pressure P and temperature conditions, $ 
is —G/T, where F is the Helmholtz energy (free energy) 
and G is the Gibbs energy (free enthalpy). These free 
entropy functions are also known as the Massieu-Planck 
potentials [12]. Usually, they are concave. 

For the finite-dimensional systems, like chemical kinet- 
ics or the Master equation, the dynamics satisfy linear 
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(linearized) kinetic equation x = Kx, where 
I 



dxidxj 



i.e. K = L(i?^$)^o 



This matrix is not symmetric but the product {D^^)x'>'iK 
= (_D^$)a;oqi(D^<i))^cq Is ahcady symmetric, hence, K is 
symmetric (self-adjoint) in the cntropic scalar product 



(a|i^6)<i, = {Ka\b)^, 



where 



(a I 6), 



dxidxj 



(2) 



(3) 



Further on, we use the angular brackets for the entropic 
inner product (j3|) and its generalizations and omit the sub- 
script $. 

For the spatially distributed systems with transport 
processes, the variables Xi are functions of the space co- 
ordinates ^, the equations of divergence form appear, 
dtXi — — • Ji, thermodynamic forces include also gradi- 
ents in space variables, Xi = V / dxi and the operator 
K has the form 



K L 

dxidxj 



i.e. K = L{D^<P)^.^A^ , 



where is the Laplace operator. This operator K is 
self-adjoint in the inner product which is just the integral 
in space of Q. The generalizations to inhomogeneous 
equilibria, non-isotropic and non-euclidian spaces are also 
routine but lead to more cumbersome formulas. 

Symmetric operators have many important properties. 
Their spectrum is real, for a function of a real variable 
/ with real values it is possible to define f{K) through 
the spectral decomposition of K, and this f{K) is also 
symmetric in the same inner product. This property is 
the cornerstone for further consideration. 

Symmetry betv^reen observables and initial data. 

— The exponential of a symmetric operator is also sym- 
metric, hence, Onsager's relations ([2]) immediately imply 



(a I eiip{Kt) b) = {b \ e^p{Kt) a) 



(4) 



The expression x{t) = exp{Kt) b gives a solution to the 
kinetic equations ^ with initial conditions x{0) = b. The 
expression {a \ x{t)) is the result of a measurement: for- 
mally, for each vector a we can introduce a "device" (an 
observer), which measures the scalar product of vector a 
on a current state x. 

The left hand side of Q represents the result of such 
an experiment: we prepare an initial state a;(0) = 6, start 
the process from this state and measure (a | x(t)). In the 
right hand side, the initial condition b and the observer a 
exchange their positions and roles: we start from the ini- 
tial condition a;(0) = a and measure (6 | x{t)). The result 
is the same function of time t. 



This exchange of the observer and the initial state trans- 
forms an ideal experiment into another ideal experiment 
(we call them dual experiments). The left and the right 
hand sides of @ represent different experimental situa- 
tions but with the same results of the measurements. 

This observation produces many consequences. As a 
first class of examples, we present the time-reversible 
Markov chains [15^, or the same class of kinetic equations, 
the monomolecular reactions with detailed balance (see 
any detailed textbook in chemical kinetics, for example, 

[S]). 

Here a terminological comment is necessary. The term 
"reversible" has three different senses in thermodynamics 
and kinetics. 

• First of all, processes with entropy growth are irre- 
versible. In this sense, all processes under considera- 
tion are irreversible. 

• Secondly, processes with microreversibility, which sat- 
isfy the detailed balance and Onsager relations, are 
time-reversible (or, for short, one often calls them 
"reversible" ) . We always call them time-reversible to 
avoid confusion. 

• In the third sense, reversibility is the existence of in- 
verse processes: if transition A ^ B exists then tran- 
sition B ^ A exists too. This condition is signifi- 
cantly weaker than microreversibility. 

"Time-reversibility" of irreversible processes sounds para- 
doxical and requires comments. The most direct interpre- 
tation of "time-reversing" is to go back in time: we take 
a solution to dynamic equations x{t) and check whether 
x{—t) is also a solution. For the microscopic dynamics (the 
Newton or Schrodinger equations) we expect that x{—t) is 
also a solution to the dynamic equations. Nonequilibrium 
statistical physics combines this idea with the description 
of macroscopic or mesoscopic kinetics by an ensemble of 
elementary processes: collisions, reactions or jumps. The 
microscopic "reversing of time" turns at this level into 
the "reversing of arrows": reaction X^i o;*^ ~^ YljPjBj 
transforms into J^jf^j-^j ~^ ^^'^ conversely. The 

equilibrium ensemble should be invariant with respect to 
this transformation. This leads us immediately to the con- 
cept of detailed balance: each process is equilibrated by its 
reverse process. "Time-reversible kinetic process" stands 
for "irreversible process with the time-reversible underly- 
ing microdynamics" . 

We consider a general network of linear reactions. This 
network is represented as a directed graph (digraph) [13] : 
vertices correspond to components Ai (i = l,2,...,n), 
edges correspond to reactions Ai — )■ Aj {i ^ j). For each 
vertex, A^, a positive real variable Ci (concentration) is de- 
fined. For each reaction, Ai Aj a nonnegative continu- 
ous bounded function, the reaction rate constant kji > 
is given. The kinetic equations have the standard Master 
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equation form 

dc 

= ^ — kjiCi) . (5) 

The principle of detailed balance ("time-reversibility") 
means that there exists such a positive vector c^'^ > 
that for all i,j [j ^ i) 

The following conditions are necessary and sufficient for 
existence of such an equilibrium cf^ > 0: 

• Reversibility (in the third sense): if kji > then kij > 
0; 

• For any cycle Ai^ Ai^ — ;> . . . -> Ai^ — > Ai^ the 
product of constants of reactions is equal to the prod- 
uct of constants of reverse reactions, 

where iq+\ — i\. This is the Wegscheider identity [5'. 

It is sufficient to consider in conditions ([7]) a finite number 
of basic cycles [T^ • 

The free entropy function for the Master equation ^ is 
the (minus) relative entropy 

In this form, the function —RTY was used already by 
L. Onsager T under the name "free energy". The entropic 
inner product for the free entropy ^ is 

(a|&)=E^- (9) 

Let c° (t) be a solution of kinetic equations ([5]) with ini- 
tial conditions c°(0) = a. Then the reciprocity relations 
^ for linear systems with detailed balance take the form 

Let us use for a and b the vectors of the standard basis 
in M": = Siq, hi = Si,., q ^ r. This choice results in 
the useful particular form of ((TU)) . We compare two ex- 
perimental situations, c°(0) = biq (the process starts at 
the state "everything is in Aq) and c^(0) = (5^^ (the pro- 
cess starts at the state "everything is in A^"); for the first 
situation we measure c"(t) ("how much Ar is produced 
from the initial Ag"), for the second one we measure Cq(i) 
("how much Aq is produced from the initial A^"). The 
reciprocal relations ([TU]) give 

eq eq • \^^) 



More examples of such relations for chemical kinetics are 
presented in |16) . It is much more straightforward to check 
experimentally these relations between kinetic curves than 
the initial Onsager relations between kinetic coefficients. 
We give an example of such an experiment below. For 
processes distributed in space, instead of concentrations 
of A and B some of their Fourier or wavelet coefficients 
appear. 

Weak form of detailed balance. For many real 
systems some of the elementary reactions are practically 
irreversible. Hence the first condition of detailed balance, 
the reversibility (if kji > then fc^ > 0) may be violated. 
Nevertheless, these systems may be considered as limits 
of systems with detailed balance when some of the con- 
stants tend to zero. For such limits, the condition ([7]) per- 
sists, and for any cycle the product of constants of direct 
reactions is equal to the product of constants of reverse 
reactions. 

This is a weak form of detailed balance without the oblig- 
atory existence of a positive equilibrium. In this section, 
we consider the systems, which satisfy this weak condition, 
the weakly time-reversible systems. 

For a linear system, the following condition is necessary 
and sufficient for its weak time-reversibility: In any cycle 
Ai-^ — >■ Ai^ —>■...—> Ai^ — > Ai-^^ with strictly positive 
constants ki^^-^i. > (here iq+i ~ ii) all the reactions are 
reversible (ki-i-^-^ > ^) <^^d, the identity ^ holds. 

The components Aq and Ar (q ^ r) are strongly con- 
nected if there exist oriented paths both from Aq to Ar 
and from Ar to Aq (each oriented edge corresponds to a 
reaction with nonzero reaction rate constant). It is con- 
venient to consider an empty path from Ai to itself as an 
oriented path. 

For strongly connected components of a weakly time- 
reversible system, all reactions in any directed path be- 
tween them are reversible. This is a structural condition 
of the weak time-reversibility. 

Under this structural condition, the classes of strongly 
connected components form a partition of the set of com- 
ponents: these classes either coincide or do not intersect 
and each component belongs to one of them. Each cycle 
belongs to one class. 

Let Aq and Aj. be strongly connected. Let us select 
an arbitrary oriented path p between Aq and Aj.: Aq 
Ai^ ^ Ai^ ^ . . . Ai^ ^ Ar. For the product of direct 
reaction rate constants in this path we use and for the 
product of reverse reaction rate constants we use . The 
ratio Krq = K^/K^ does not depend on the path p and 
characterizes the pair Ar,Aq, because of the Wegscheider 
identity This is the quantitative criterion of the weak 
time-reversibility. 

The constant Krq is an analogue to the equilibrium con- 
stant. Indeed, for the systems with positive equilibrium 
and detailed balance, KrqC^^ — c^^ and Krq — c^'^/c"'^. 

For weakly time-reversible system, the reciprocal rela- 
tions between kinetic curves can be formulated for any 
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strongly connected pair Aq and Ar. Exactly for the same 
pair of kinetic curves, as in (|lip. we obtain 



The analytic solution easily gives 



cm 



Kr 



(12) 



This formula describes two experiments: (i) we start the 
system at t = from the pure Aq and measure Cr{t), then 
(ii) we start at t = from the pure Ar and measure Cq{t). 
The ratio of these two kinetic curves, Cr{t)/cq{t) does not 
depend on t and is equal to the generalized equilibrium 
constant K^-q- 

The weak form of the Wegscheider identity for general 
(nonlinear) kinetic systems is also possible. Let us con- 
sider the reaction system: 

arlAi + . . . + arnAn PrlAi + . . . + /3rn^n , (13) 

which satisfies the mass action law: c = ^^kr Y\i c"' j 
where kr > 0, ^ri — Pri ~ ctri is the stoichiometric vector 
of the rth reaction, and the reverse reactions with positive 
constants are included in the list ([T^ separately. 

Let us consider linear relations between vectors {jr}- 



Ar7r = and 7^ for some r . 



(14) 



If all the reactions are reversible then the principle of de- 
tailed balance gives us the identity |T3]: 



(15) 



for any linear relation (|14p . For reversible reactions, we 
can take Ar > in (IT5t for all r: if we substitute the 
reactions with A^ < by their reverse reactions then 7^ 
and Ar change signs. It is sufficient to consider only the 
cone A+ of non-negative relations (|14p (A^ > 0) and take in 
(fT5)) the direction vectors of its extreme rays. Let = 
for some r. The weak form of the identity dTSl) is: 

For any extreme ray 0/ A+ with a direction vector A^ > 
the reactions which correspond to the positive coejficients 
Ar > are reversible (k~ > Q) and their constants satisfy 
the identity hl5\l . 

Nonlinear Examples. — It seems impossible to find 
a general relation between kinetic curves for general non- 
linear kinetics far from equilibrium. Nevertheless, sim- 
ple examples encourage us to look for a nontrivial theory 
for some classes of nonlinear systems. In this Section, we 
give two examples of nonlinear elementary reactions which 
demonstrate the equilibrium relations between nonequilib- 
rium kinetic curves |16j . 

2 A o B. The linear conservation law is Cyi -|- 2cs = 
const. Let us take two initial states with the same value 
CA + 2cB - 1: (a) ca(0) = 1, cs(0) = and {b) ca(0) = 
0, cb(0) = 1/2. We will mark the corresponding solutions 
by the upper indexes a,b. The mass action law gives: 

CA - -2fc+4 + fc-(l - ca) , cb = (1 - ca)I2 . (16) 



k+_ 
IF 



(17) 



the denominator involves the A concentrations of both 
trajectories, c° (started from c^(0) = 1, cb(0) = 0) and c^ 
(started from ca{0) = 0, cb{0) — 1/2). A ratio is equal to 
the equilibrium constant at every time t > 0. This identity 
between the non-stationary kinetic curves reproduces the 
equilibrium ratio. 

2 A o 2B. The linear conservation law is ca + cb = 
const. Let us take two initial states with the same value 
CA + CB ^ 1: (a) ca(0) = 1, cb(0) = and (&) ca(0) = 
0, cb{0) — 1. The kinetic equation is 



CA 



-2k+c\ + k-{l-CAf , CB = l-CA. (18) 



It can be solved analytically. For this solution. 



C|(^)4(t) ^k^^ ^^.eq ^ 

c^(t)4(t) k- 



(19) 



both the numerator and denominator include trajectories 
for both initial states, a and b. This identity between the 
kinetic curves also reproduces the equilibrium ratio. 

Experimental evidences. — In this work, we investi- 
gate the validity of the reciprocal relations using the TAP 
(Temporal Analysis of Products) technique proposed by 
Cleaves in 1988 jTT]. It has been successfully applied 
in many areas of chemical kinetics and engineering for 
non-steady-state kinetic characterization [H]. The stud- 
ied reaction is a part of the reversible water gas shift re- 
action over iron oxide catalyst. The overall reaction is 
H2O + G0^ H2 + CO2. 

Experimental set-up. The TAP reactor system used 
in this work is made of quartz and is of the size 33 mm 
bed-length and 4.75 mm inner diameter. The products 
and the unreacted reactants coming out of the reactor are 
monitored by a UTI lOOC quadrupole mass spectrometer 
(QMS). The number of molecules admitted during pulse 
experiments amounts to 10"'^^ molecules/pulse. 

To ensure uniformity of the catalyst along the bed, we 
use a thin-zone TAP reactor (TZTR), the width of the 
catalyst zone being 2mm. Experiments were performed 
over 40 mg of Fe203 catalyst. The catalyst was packed 
in between two inert zones of quartz particles of the same 
size (250 < dp < 500/xm). The temperature of the reactor 
was measured by a thermocouple positioned in the center 
of the catalyst bed. Several single pulse experiments were 
performed by pulsing CO or CO2 at the temperature of 
780K. In all the experiments, the reaction mixture was 
prepared with Ar as one of the components, so that the 
inlet amount of the components can be determined from 
the Ar response. 
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Fig. 1: Fourier domain result values for the "B from A/ A from 
B" ratio H20|) . vs. frequency / in Hz (so that lj = 2nf, s = iio); 
real and imaginary part. The error bars were obtained from 
10,000 resampled measurements. 

Application to the measurements. In a thin-zone 
TAP-reactor, the diffusion occurring in the inert zones 
flanking the thin reactive zone must be accounted for. 
The Knudsen regime in these zones guarantees a linear 
behaviour, so that the resulting outlet fluxes can be ex- 
pressed in terms of convolutions. Switching to the Laplace 
domain greatly facilitates the analysis, and we can prove 
in general that the fixed proportion property is equivalent 
to the following equality in terms of the exit fluxes -Fba 
of gas B given a unit inlet pulse of gas A and , of A 
given a unit inlet pulse of B, see |19) : 

^ ^ (cosh y^sn;X){y^ sinh y^sTXB) CFb^ js) ^2^^ 
(cosh ysTTs ) (y/TXA sinh ^/stxa) CFab («) 

identically in the Laplace variable s, where t^^g = 
CiL^/Do, Ci denoting the packing density of the i-th zone, 
Li its length, and Do the diffusivity of gas G. To apply 
this in practice, we set s = iut and switch to the Fourier 
domain. 

Performing these corrections, with A denoting CO and 
B CO 2, the results of Fig. [T] in the Fourier domain are 
obtained. The real and imaginary parts of the right-hand 
side in (|20p are graphed, with error bars corresponding to 
three times the standard deviation estimated from resam- 
pling 10,000 times the exit flux measurements using their 
principal error components. Ideally, all imaginary values 
should be zero; we see that zero does lie within all the 
confidence intervals. We also see that the smallest error 
in the real parts occurs for the second frequency, 2.2 Hz. 
This confidence interval lies snugly within the others, of- 
fering confirmation that (within experimental error) the 
same value for all frequencies is obtained. 

Conclusion. — The shift in time operator is symmet- 
ric in the entropic inner product. Its symmetry allows us 
to formulate the symmetry relations between the observ- 
ables and initial data. These relations could be validated 
without differentiation of empiric curves and are, in that 



sense, more robust and closer to the direct measurements. 
For the Markov processes and chemical kinetics, the sym- 
metry relations between the observables and initial data 
have an elegant form of the symmetry between "A pro- 
duced from i?" and "i? produced from A": their ratio is 
equal to the equilibrium constant and does not change in 
time pT|). For processes distributed in space, in- 

stead of concentrations of A and B some of their Fourier 
or wavelet coefficients appear. 

The symmetry relations between the observables and 
initial data have a rich variety of realizations, which makes 
the direct experimental verification possible. On the other 
hand, this symmetry provides the possibility to extract 
information about the experimental data through the dual 
experiments. These relations are applicable to all systems 
with microreversibility. 
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